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Mirror dark matter interacting with ordinary matter via photon-mirror photon kinetic 
mixing can explain the DAMA, CoGeNT and CRESST-II direct detection experiments. 
This explanation requires kinetic mixing of strength e ~ 10 -9 . Such kinetic mixing will 
have important implications for early Universe cosmology. We calculate the additional 
relativistic energy density at recombination, SN e ff[CMB]. We also calculate the effects 
for big bang nucleosynthesis, 5N e ff[BBN]. Current hints that both 5N e ff[CMB] and 
SN e ff[BBN] are non-zero and positive can be accommodated within this framework if 
e ~ few x 1CT 9 . In the near future, measurements from the Planck mission will either 
confirm these hints or constrain e ~ 10 -9 . 
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The long standing DAMA annual modulation signalfl], along with more recent data 
from the CoGeNT[2] and CRESST-II[3] experiments can be explained within a simple 
hidden sector framework[4, 5] (see also ref.[6, 7] for earlier studies). The requirement 
is for a hidden sector with an unbroken U(l)' interaction kinetically mixed with the 
standard U(l)y gauge field. Such a framework can explain all of the direct detection 
results, both positive and negative, provided that the hidden sector contains two or more 
stable particles. Mirror dark matter [8] appears to be the most well motivated and also 
constrained example of such a theory (see ref.[9] for reviews and more comprehensive 
references). In the mirror dark matter theory, the hidden sector is exactly isomorphic to 
the ordinary sector except the chiralities are interchanged. In this case an exact unbroken 
mirror symmetry can be defined which interchanges each ordinary particle with a mirror 
partner, as well as mapping x, t — > —x, t[8}. The latter property allows mirror symmetry to 
be identified with space-time parity which provides an interesting theoretical motivation 
for such a theory 2 . 

In this theory dark matter consists of a spectrum of stable mirror particles of known 
masses (e', H', He', O', ...) which interact with the ordinary particles via gravity and 
renormalizable photon-mirror photon kinetic mixing: [12] 3 

'-mix 2 fxv K 1 ) 

where F^ is the ordinary photon field strength tensor, and F' is the field strength tensor 
for the mirror photon. This interaction enables mirror charged particles, such as e', p', to 
couple to ordinary photons with electric charge ee [14]. 

Studies have shown [15] that mirror dark matter can provide a successful account of 
the Large Scale Structure (LSS) of the Universe which can mimic cold dark matter or 
warm dark matter, depending on initial conditions. Mirror dark matter, though, has 
quite distinctive properties on smaller scales. Observations suggest that a) In galaxy 
clusters or at least in some of them, mirror dark matter should be predominately confined 
within galactic halos to account for the observations of the Bullet cluster[16] (c.f. [17]). 
b) Within galactic halos, limits on the MACHO fraction [18] suggest that mirror dark 
matter is mainly (at least ~ 50%) in gaseous form - i.e. an ionized plasma comprising e', 
H', He', O', .... A substantial MACHO fraction is also possible[19]. Mirror dark matter 
is dissipative, which indicates that a substantial heat source should exist to replace the 
energy lost due to radiative cooling. It turns out that ordinary supernova can supply 
the required energy to stabilize the mirror dark matter halo provided that e ~ 1CT 9 [20]. 
This heat source is anisotropic in spiral galaxies (given that the supernova originate 
predominately from the disk) which might explain the deviations from perfect spherical 
halos which are necessary to agree with observations, c) The self-interactions of mirror 

2 There are also interesting models with mirror symmetry spontaneously broken, see e.g. [10, 11]. 
3 Other interactions are possible including a Higgs-mirror Higgs coupling Xcj)^ (jxj)'^ <j)' , which can lead to 
novel effects for Higgs physics that can be probed at the LHC[13]. 



1 



dark matter within galactic halos might help explain the inferred cored dark matter 
distribution - a long standing puzzle of the standard cold dark matter paradigm. 

The gaseous component of the halo is presumed to be composed of an ionized plasma 
comprising a spectrum of components e', H', He', O', These components can po- 
tentially be detected in direct detection experiments via the kinetic mixing interaction, 
Eq.(l), which induces spin-independent elastic Rutherford scattering. The H', He' com- 
ponents, which are naively expected to be dominant, are however very light. In fact 
these components are so light that they are unable to explain the DAMA, CoGeNT or 
CRESST-II signals. The experiment most sensitive to the He' component turns out to be 
the earlier CRESST/Sapphire experiment [21]. That experiment can be used to infer the 
bound: [6] 



5 3 x 1(T 9 . (2) 



The on-going Texono experiment [22] at Jin-Ping underground laboratory is expected to 
provide a much more sensitive probe of the He' component in the near future. The 
remaining heavier components, are for simplicity assumed to be composed of a 'dominant' 
metal component, A'. The direct detection experiments, DAMA, CoGeNT and CRESST- 
II can be explained by the A' component if [4, 5]: 



e\J^A' ~ few x 10 
m A > 



-10 



16 - 56 (3) 



where £4/ = n^/m^/ / '(0.3 GeV/cm 3 ) is the halo mass fraction of the species A' and m p 
is the proton mass. Combining Eq.(2) and Eq.(3), we obtain a lower limit on the ratio 
^A'/^He' ~ 10~ 2 at the Earth's location in the halo. This appears to be a relatively 
high metal component, and might be due to rapid mirror star formation and evolution 
at an early epoch [23]. It might also be possible to produce the metal component in the 
early Universe, although this is disfavoured in the simplest scenarios with high reheating 
temperature [24]. Allowing for possible uncertainties in £4/, Eq.(3) suggests that e should 
be in the range: 

10- 10 £ e £ 10- 8 . (4) 

This range of e is consistent with the direct laboratory limits, the most stringent of which 
arises from invisible decays of orthopositronium[25, 26], e ~ 1.5 x 10~ 7 . An important 
proposal exists [27] for a more sensitive orthopositronium experiment which can cover 
much of the e range of interest, Eq.(4). 4 

In the context of early Universe cosmology, the kinetic mixing can induce processes 
which transfer energy into the mirror sector. This will lead to various modifications to 



4 There are a number of other possible implications of kinetic mixing, for example if the solar system 
contains mirror matter space-bodies [28], see also the review [29]. 
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big bang nucleosynthesis (BBN), Cosmic Microwave Background anisotropies (CMB) and 
also Large Scale Structure (LSS). This allows values of e ~ 1CT 9 to be probed in early- 
Universe cosmology. 

Recall that the relativistic energy component at recombination can be parameterized 
in terms of the effective number of neutrino species, N e ff[CMB], by: 



where p 7 is the energy density of the CMB photons. In the standard electroweak theory, 
N e ff ~ 3.046 (the difference from 3 being due to slight heating of the neutrinos from 
ee annihilation) [30] . Currently there are some interesting hints that 5N e ff[CMB] = 
N e ff[CMB] — 3.046 is in fact non-zero and positive within the usual ACDM frame- 
work. For example, the Wilkinson Microwave Anisotropy Probe (WMAP) obtained 
N eff [CMB] = 4.34 ± 0.87 (1 a C.L.) [31]. Observations with the Atacama Cosmol- 
ogy Telescope suggest N eff [CMB] = 4.56 ±0.75 (1 a C.L) [32]. Also, results of the South 
Pole Telescope combined with WMAP7 suggest N eff [CMB] = 3.86 ±0.42 (1 a C.L.) [33]. 
A recent analysis [34] combining the data finds that: 



An effective neutrino number can also be defined for BBN, and it is worth noting that 
two recent measurements [35, 36] of Y p suggest N e ff[BBN] > 3. These results have moti- 
vated a number of proposals for new physics, including models with sterile neutrinos [3 7]. 
Importantly, high precision results from the Planck mission should be able to confirm a 
non-zero 5N e ff[CMB] to within a precision of around 0.2[38] in the near future. In the 
meantime, it is timely to work out the predictions for N e ff[CMB] and N e ff[BBN] in the 
mirror dark matter model. 

In the e — > limit, the ordinary and mirror particles decouple from each other, and 
in general they may have different temperatures, T, T'. It will be useful to define the 
temperatures: T 7 [T 7 ] for the temperature of the ordinary [mirror] photons and T u \T' U ] 
is the temperature of the ordinary [mirror] neutrinos. In this study, we assume the 
initial conditions T' ,T' V <C T 7 = T v due to some physics at early times. Asymmetric 
reheating within inflationary scenarios is one possibility [39]. If we assume that only the 
ordinary matter is reheated after inflation, then we have the initial condition T' ,T' V ~ 
and T 7 = T v = T RH - With these initial conditions it is always safe to ignore T' v since 
T' u <C T' f in the period of interest. This is because T 7 /T 7 generation occurs mainly in the 

low temperature region, T ~ few MeV, where the mirror weak interaction rate is always 
much less than the expansion rate: G 2 F T' 5 <C T 2 jm p \. 

For T 7 < 100 MeV, ee — > eV is the dominant process which transfers energy into the 
mirror sector. This process will not only generate but will also slightly reduce T 7 . This 
energy transfer happens predominately after the kinetic decoupling of the neutrinos (i.e. 





N eff [CMB] = 4.08 



+0.71 

-0.68 



at 95% C.L. 



(6) 
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for T 7 ~ 3 MeV), so the effect is to slightly decrease T y /T v . Therefore, at recombination 
the energy density of neutrinos will be larger than it would otherwise have been if there 
were no kinetic mixing. This effect will contribute to SN e ff[CMB] along with the more 
obvious contribution to 5N e ff[CMB] from the increase in energy density due to the 
production of the mirror photons. We will see that these two contributions are numerically 
similar in magnitude. 

The cross-section for the process ee — > e'e' is: 

a=^a^- 3 (s + 2mlY (7) 

o S 

where s is the usual Mandelstam variable and m e is the electron mass. The energy transfer 
to the mirror sector within a co-moving volume, R 3 , per unit time is: 

^ = R 3 n e +n e -(av M 0i£) (8) 

where 8 is the energy transferred in the process ee — > e'e', v m0 i is the M0ller velocity (see 
e.g. ref.[40]), and n e - ~ n e + is the number density of electrons: 



^/E 2 - m\E 
+ exp{E/T 1 ) 

The quantity {av M 0i£} is given by [41] (see also ref.[40]): 



1 r oo JEP-mlE 
= J_ / _V dE . ( 9 ) 

7T 2 J me l + exv(E T~) v 1 



oo f 00 f°° / E 2 

{ ^- £} = ^ptff L ds ° (s - im '^L ^- e ^m^ - 1 

(10) 

where K 2 is the modified Bessel function of order 2. The quantity oo 0.8 takes into 
account various approximations in deriving Eq.(10) such as replacing the exact Fermi- 
Dirac distribution of e, e with the simpler Maxwellian one[41]. 

Energy is transferred from the ordinary particles to the mirror particles. The second 
law of thermodynamics can be used to work out the change in entropy of the ordinary 
particles: 

iS = ^. (11) 

To a very good approximation, we can neglect the contribution of the neutrinos to dS 
since energy is mostly transferred to the mirror sector after neutrino kinetic decoupling 
(i.e. for T 7 ~ 3 MeV). Thus, R oc 1/T V , and the above equation has the form: 



d_ 

dt 



(p 7 +P 7 + p e +Pe)R 3 



T 

7 



n e -n e + (av mo i£)R 3 
T 7 



4 



where 



Pi 

Pe 
Pe 



15 7 
Pi 

3 

7T 2 

3tt 2 



( u 2 _ ^2)1/2^2 



1 + e u 



du 



{u 2 



X 



2N3/2 



1 + e u 



du 



and x = m e /T 7 . The entropy gained by the mirror particles is: 



V 

7 



This equation implies that 



eft 



(P 7 +P 7 + Pe+Pe)^ 3 

T' 

7 



n e -n e +(<Jv mol £)R 3 
V 

7 



where 



^7 = 



^ T ,4 

15 7 
3 



Pe 



2T^ roo ( u 2_ a ,/2 )1/2w2 



_ zi 7 r 

7T 2 is' 



1 + e" 



2T: 



7-1 



3tt 2 



it 



/2N3/2 



1 + e u 
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and x' — m e /T'. The Friedmann equation reads: 



(13) 



(14) 



(15) 



(16) 



R 
R 



8tt 
3m 2 ; 



Pj + Pe + Pv + P., + P, 



(17) 



where m p i ~ 1.22 x 10 22 MeV is the Planck mass. The energy density of the neutrinos is 
given by the standard value, p u = Eq.(12,15) and Eq.(17) form a closed system 

which can be solved to give the evolution of T 7 , T u and T' . By way of example, the result 
for e = 10 -9 is shown in figure 1. 
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Figure 1: Evolution of T 7 (solid line), T u (dashed line) and T' (dotted line) for e = 10 9 . 



We find numerically that T'/T 7 evolves to a constant which satisfies: 



V 

-1 -0.31 
T 



10- 9 



1/2 



(18) 



which is consistent with the results of the earlier study [41]. 

The energy density of the relativistic components at recombination can be expressed 
in terms of an effective neutrino number, Eq.(5). The change in this effective neutrino 
number due to the kinetic mixing effect, 5N e ff[CMB], can be computed as follows: 



6N eff [CMB] = 3 



T 



T u (e = 0) 



- 1 + 



7 V71,(e = 0), 



(19) 



where the temperatures are evaluated at the time when photon decoupling occurs, i.e. 
when T 7 = Td ec = 0.26 eV. The first term on the right-hand side, 8N£ff, is the contribution 
due to the photon cooling, which increases the relative energy density of the neutrinos. 
The second term, dN^f, is the contribution to the energy density from the mirror photons. 
In figure 2 we plot 5N e ff[CMB] versus e, showing also the separate contributions 8N£fj 
and SN^ff. As remarked earlier [see the discussion around Eq.(6)] there are current hints 
that SN e ff[CMB] is non-zero and positive. Indeed, the result of the combined analysis[34], 
Eq.(6), provides an interesting hint that e few x 10~ 9 . 
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Figure 2: 6N e ff[CMB] versus e (solid line). The separate contributions, discussed in the text, 
<W°// (dashed line) and 5N^^ (dotted line) are also shown. 

Previous studies[15] based on successful LSS have suggested that T' /T^ should obey an 
upper limit of around 0.3. As pointed out previously [41], Eq.(18) then suggests an upper 
limit on e of around e ~ 1CT 9 . However due to the complexity of modeling the effects of 
mirror dark matter for LSS there is a significant uncertainty (maybe even a factor of two 
or so) in the limit on T 7 /T 7 from the existing studies based on LSS considerations, and 

thus a more conservative upper limit on e is likely to be around e ~ 4 x 1CT 9 . Obviously 
if the Planck mission confirms the current hints that SN e ff[CMB] > 0, then this would 
warrant a careful re-analysis of the LSS implications of the mirror dark matter theory. 

The transfer of energy into the mirror sector also affects BBN as first pointed out 
in ref. [42]. The effect on the primordial Helium abundance can easily be obtained by 
computing Y p in the usual way, but with the T 7 , T v dependent rates: n + e <H> p + P e , n + 
v e ■<->■ p + e, n ■<-)■ p + e + P e evolved down to the deuterium 'bottle neck' temperature 
T 7 = 0.07 MeV obtained by utilizing the solution of Eq. (12,15) and Eq.(17) 5 . In this way 
we can compute the helium mass fraction for a particular value of e, Y p (e). For e = 0, the 

5 BBN implications of mirror dark matter without kinetic mixing, but with non-zero T 1 /T arising from 
assumed initial conditions was studied in ref. [43]. 
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result is the standard value, Y p (0) ~ 0.24. We define SN e ff[BBN] in the usual way: 

5N eff [BBN) = Ypie) 0M l M • (20) 

Thus 5N e ff[BBN] = 1 corresponds to an increase of Y p of 0.013 which is equivalent to 
having 4 neutrinos instead of 3 contributing to the energy density. In figure 3 we plot 
our results for 5N e ff[BBN] versus e. Comparison of figure 3 with figure 2 shows that we 
expect 5N e ff[CMB] > SN e ff[BBN}. The diminished effect for BBN is simply because 
a significant part of the generation of and reduction of T 7 /T^ occurs below the freeze 
out temperature of weak interactions T w ps 0.8 MeV. 

1.2 i 1 1 1 1 1 1 1 

1 - 



0.8 - 
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£ 

Figure 3: 5N eff [BBN] versus e. 

In conclusion, we have computed the additional radiation content of the Universe 
at the recombination epoch due to the effects of photon-mirror photon kinetic mix- 
ing, 5N e ff[CMB]. We have also calculated the effects for big bang nucleosynthesis, 
5N eff [BBN}. There are some interesting hints that both 5N eff [CMB] and 5N eff [BBN] 
are non-zero and positive, which if confirmed by e.g. the anticipated high precision mea- 
surements of the CMB from the Planck mission, will ultimately allow e to be determined. 
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